A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

An identification method to estimate surface shortwave fluxes on temperature and thickness observations of sea ice in CHINARE 2006

In this paper, a parameter identification method to determine surface shortwave fluxes using temperature and thickness measurements of sea ice in CHINARE 2006 is presented. Adopting a new standard for the calculation of the thermodynamic properties of seawater named TEOS‐10, the surface shortwave fluxes are calculated by the temperature and thickness observations that were measured at Nella Fjord around Zhongshan Station, Antarctica. New simulations for temperature profiles in a different measurement period are performed by three parameterization schemes including the present method, Zillman and Shine. All numerical results are compared with in situ measurements. Results show that better simulations of the surface shortwave radiations and temperature distributions are possible with the identification method than Zillman and Shine. Therefore this method is valid, and the obtained shortwave radiation function can be applied in sea ice modeling. Copyright © 2015 John Wiley & Sons, Ltd.     

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

一种提高极值点处收敛精度的三阶WENO-Z格式

为了提高三阶WENO-Z格式在极值点处的计算精度，通过理论推导给出三阶WENO格式满足收敛精度的充分条件。采用泰勒级数展开的方式，推导给出所构造格式非线性权重的计算公式，并综合权衡计算精度和计算稳定性确定所构造格式的参数。通过两个典型的精度测试，验证了改进格式在光滑流场极值点区域逼近三阶精度。进一步选用激波与熵波相互作用和Richtmyer-Meshkov不稳定性等经典算例，证实了本文提出的改进格式WENO-PZ3相较其他格式（WENO-JS3和WENO-Z3）不仅具有较高的精度，而且降低了格式的耗散，提高了对流场结构的分辨率。     

Water waves trapped by thin submerged cylinders in a two‐layer fluid: Discrete eigenvalue

Exact solutions of the linear water‐wave problem describing oblique water waves trapped by a submerged horizontal cylinder of small (but otherwise fairly arbitrary) cross‐section in a two‐layer fluid are constructed in the form of convergent series in powers of the small parameter characterising the “thinness” of the cylinder. The terms of this series are expressed through the solutions of the exterior Neumann problem for the Laplace equation describing the flow of unbounded fluid past the cylinder.     

The impact of a conceptual model-based mathematics computer tutor on multiplicative reasoning and problem-solving of students with learning disabilities

The study explored the impact of Please Go Bring Me-COnceptual Model-based Problem Solving (PGBM-COMPS) computer tutoring system on multiplicative reasoning and problem solving of students with learning disabilities. The PGBM-COMPS program focused on enhancing the multiplicative reasoning and problem solving through nurturing fundamental mathematical ideas and moving students above and beyond the concrete level of operation. This is achieved by taking advantages of the constructivist approach from mathematics education and explicit conceptual model-based problem solving approach from special education. Participants were three elementary students with learning disabilities (LD). A mixed method design was employed to investigate the effect of the PGBM-COMPS program on enhancing students’ multiplicative reasoning and problem solving. It was found that the PGBM-COMPS program significantly improved participating students’ problem solving performance not only on researcher developed criterion tests but also on a norm-referenced standardized test. Qualitative and quantities data from this study indicate that, in addition to nurturing fundamental concept of composite units, it is necessary to help students to understand underlying problem structures and move toward mathematical model-based problem representation and solving for generalized problem solving skills.     

CH4/空气简化动力学机理数值验证

选择绕圆柱预混燃烧算例，验证CH₄/空气三种简化动力学机理（16s41r、15s19r和53s325r）.考虑均匀来流，忽略湍流和湍流与燃烧相互作用以及燃料扩散效应，假设层流有限反应速率，采用保自由流5阶WENO格式求解多组分Euler方程组，得到CH₄/空气预混燃烧流场温度等值线、沿驻点线压力和温度及其CH₄、CO和CO₂质量百分数分布.结果表明：三种简化动力学机理给出的流场均出现弓形激波和火焰面，弓形激波和火焰驻点距离及其形状、诱导区宽度和简化动力学机理相关.当圆柱直径增大，弓形激波和火焰向圆柱上游移动，对应的驻点距离均增大，诱导区宽度变短，点火延时变小，但火焰和弓形激波位置次序未变化.53s325r模型要比16s41r模型和15s19r模型精度要高，点火延时覆盖的压力和温度范围也较宽，所有简化机理均未完全反应，在较大圆柱直径下游达到化学平衡.     

Hyperplane sections of determinantal varieties over finite fields and linear codes

We determine the number of ${\mathbb{F}}_q$-rational points of hyperplane sections of classical determinantal varieties defined by the vanishing of minors of a fixed size of a generic matrix, and identify the hyperplane sections giving the maximum number of ${\mathbb{F}}_q$-rational points. Further we consider similar questions for sections by linear subvarieties of a fixed codimension in the ambient projective space. This is closely related to the study of linear codes associated to determinantal varieties, and the determination of their weight distribution, minimum distance, and generalized Hamming weights. The previously known results about these are generalized and expanded significantly. Connections to eigenvalues of certain association schemes, distance regular graphs, and rank metric codes are also indicated.     

Solutions with concentration for conservation laws with discontinuous flux and its applications to numerical schemes for hyperbolic systems

Measure-valued weak solutions for conservation laws with discontinuous flux are proposed and explicit formulae have been derived. We propose convergent discontinuous flux-based numerical schemes for the class of hyperbolic systems that admit nonclassical -shocks, by extending the theory of discontinuous flux for nonlinear conservation laws to scalar transport equation with a discontinuous coefficient. The article also discusses the concentration phenomenon of solutions along the line of discontinuity, for scalar transport equations with a discontinuous coefficient. The existence of the solutions for transport equation is shown using the vanishing viscosity approach and the asymptotic behavior of the solutions is also established. The performance of the numerical schemes for both scalar conservation laws and systems to capture the -shocks effectively is displayed through various numerical experiments.     

Development of LBGK and incompressible LBGK‐based lattice Boltzmann flux solvers for simulation of incompressible flows

This paper presents lattice Boltzmann Bhatnagar–Gross–Krook (LBGK) model and incompressible LBGK model‐based lattice Boltzmann flux solvers (LBFS) for simulation of incompressible flows. LBFS applies the finite volume method to directly discretize the governing differential equations recovered by lattice Boltzmann equations. The fluxes of LBFS at each cell interface are evaluated by local reconstruction of lattice Boltzmann solution. Because LBFS is applied locally at each cell interface independently, it removes the major drawbacks of conventional lattice Boltzmann method such as lattice uniformity, coupling between mesh spacing, and time interval. With LBGK and incompressible LBGK models, LBFS are examined by simulating decaying vortex flow, polar cavity flow, plane Poiseuille flow, Womersley flow, and double shear flows. The obtained numerical results show that both the LBGK and incompressible LBGK‐based LBFS have the second order of accuracy and high computational efficiency on nonuniform grids. Furthermore, LBFS with both LBGK models are also stable for the double shear flows at a high Reynolds number of 10⁵. However, for the pressure‐driven plane Poiseuille flow, when the pressure gradient is increased, the relative error associated with LBGK model grows faster than that associated with incompressible LBGK model. It seems that the incompressible LBGK‐based LBFS is more suitable for simulating incompressible flows with large pressure gradients. Copyright © 2014 John Wiley & Sons, Ltd.